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20.  (contd) 

latitudes  |A|  <  65°  (that  is,  for  Mcilwain  parameters  L  <  5.  6  re)  the  cutoff 
rigidities  are  most  sensitive  to  variations  of  the  z-component  of  the  magnetic 
field  within  geocentric  distances  r  <  2  I.recos2»  and  within  a  longitude  sector 
extending  from  the  longitude  of  the  station  up  to  60°  to  90°  to  the  east  of  it. 
Within  this  region,  magnetic  perturbations  being  overhead  and  immediately 
to  the  east  of  the  station  are  found  to  have  the  largest  relative  effect. 

In  the  second  part  of  the  study,  the  effects  of  magnetosphe ric  current 
systems,  corresponding  to  the  Bostrom  models,  on  the  cutoff  rigidities  and 
the  asymptotic  directions  at  Kiel,  .Jungfraujoch  and  Home  were  examined. 

The  effect  of  a  Bostrom  type  1  model  current  system  (partial  ring  current)  on 
cosmic  ray  cutoff  rigidities  and  asymptotic  directions  is  largest  at  longitudes 
corresponding  to  the  western  half  of  the  region  covering  the  main  disturbances 
of  the  equatorial  surface  magnetic  field.  The  dependence  on  longitude  is  more 
pronounced  at  high  latitudes  than  at  mid-  and  low-latitudes.  At  latitudes 
1  A |  >  52°  the  magnitude  of  the  effect  depends  significantly  on  the  radius  of 
the  partial  ring  current.  At  rigidities  a  few  GV  above  the  main  cutoff  the 
effect  on  asymptotic  directions  is  mainly  restricted  to  a  westward  stiift  in 
longitude  with  amplitudes  <25°  for  a  current  intensity  corresponding  to 

a  strong  magnetic  storm.  At  lower  rigidities  the  magnitude  of  the  effect 
increases  and  the  asymptotic  latitudes  are  substantially  affected  us  well.  As 

1  an  interesting  fact,  however,  the  asymptotic  direction  corresponding  to  the 
main  cutoff  rigidity  is  found  to  be  almost  invariant,  even  if  the  main  cutoff 
itself  may  have  changed  significantly.  Bostrom  No.  2  model  currents  can 
also  considerably  affect  the  cosmic  ray  cutoff  rigidities  as<  well  as  the 
asymptotic  directions.  The  maximum  effect  occurs  immediately  outside  the 
eastern  and  western  border  /ones.  The  maximum  amplitude  of  the  effect 
decreases  approximately  as  an  exponential  function  of  the  latitudinal  distance 
between  the  point  of  observation  and  the  foot  point  of  the  current  system  in  the 
ionosphere. 

The  significance  of  the  results  is  discussed  with  respect  to  cosmic  ray 
research  as  well  as  to  the  study  of  the  perturbed  magnetic  field  in  the  | 

magnetosphe  re. 
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Effects  of  Asymmetric  Magnetospheric 
Currents  on  Cosmic  Radiation 


I.  IM'KdlHCTION 

During  1  sj. •  ■  ’i : i  •  gin storms  significant  incrt'ascs  in  the  ground-level  cos¬ 
mic  ni  mti.Misit'.  otic  to  ; t  depression  of  the  cosmic  ray  cutoff  rigidities  have  been 

observed  with  .  n  adit  tides  depending  both  on  geomagnetic  latitude  and  local  time 
1  ’’ 

(Arens;  Debrunner  el  1;“  and  references  therein).  It  is  assumed  that  this  effect 
originates  from  rurrent  s\  steins  developing  in  the  earth's  magnetosphere  during 
tlie  main  phase  of  .  geom  .gnrtic  storm.  Accordingly,  several  attempts  have  been 
made  to  correlate  the  variations  of  the  eosmie  ray  cutoff  rigidities  with  tlu'  param¬ 
eters  of  v  1 1  tons  current  S'  stems.  •  I  p  to  now ,  most  of  these  studies  w  e  re 

■  ■irrjed  out  I.-.. nmitig  a  a  in  rit  ha.  11s  s\mmetric  eurrents,  neglect  ing  the  well  estab¬ 
lished  longitudinal  isv  mmet  rs  in  *!ie  uromngnet  iv  disturham-e  at  low-  and  mid  -  lat  it  tides 
dn  tint;  in  tgnefic  storm...  1  lie  ■  a  mine',  rie  perturbation  superposed  upon  a  sy  m  mot  lie 
field  depre  s  s  j.  m  1 1  a  turn  •  1 1  rib' i*  •  I  la,  si  \  c  r  d  a  ut  hors  to  i  long  it  ud  ina  11  v  limited, 

1 1  '  rt  i  !  i;nr  >  a  rent  it !  i  field  ■■ !  n  ue.  I  .  1  •  •  —  i  "e  1 1 1  roi  igh  the  mi'or  i!  i  nuns  pile  re 
f»  , 

! (  nr  riing.s.  i  ni;::  -  him  •  and  I.  m  i  I.  .  '  a  1  •  :  no  ..  ,  , i.  |t  h  ,s  been  «he 

'  i  report  *o  .— 1 1 1  <  I  1 1 1 1  ■  offer's  id'  ,-seii  m  ig- 


nlllclia 


■cl  |a. 


domain  of  the  earth’s  magnetosphere.  The  corresponding  results  were  expected 
to  demonstrate  the  significance  of  cosmic  ray  measurements  as  an  additional  tool 
for  magnetospheric  studies,  especially  in  modeling  the  disturbed  magnetic  field 
in  the  magnetosphere. 

in  the  following,  first  the  approach  to  the  problem  is  reviewed  in  Section  2, 
and  examples  of  the  basic  tools  needed  in  the  procedure  are  illustrated  in  Section  3. 
Then,  in  Section  4,  the  results  of  a  simple  study  correlating  local  perturbations 
in  the  geomagnetic  field  with  cosmic  ray  cutoff  rigidity  variations  are  discussed. 
Section  5  contains  a  summary  of  a  detailed  quantitative  analysis  on  the  effects 
of  a  partial  ring  current  and  of  field-aligned  currents  on  both  cutoff  rigidities  and 
asymptotic  directions  in  different  latitude  regions.  In  conclusion,  the  significance 
of  the  results  obtained  in  the  entire  study  is  discussed  in  Section  6  for  both  mag- 
netospheric  and  cosmic  ray  research. 


2.  METHOD  OP  PROCEDURE 


The  differential  equation  describing  the  path  of  a  particle  of  charge  q  and 
momentum  p  in  the  earth’s  magnetic  field,  li,  is  given  by 


■>-<  -» 


(1) 


(Rossi  andOlbort^),  where  x  represents  the  position  vector,  s  the  arc  length  along 
the  trajectory,  andp  =  |p|  .  Unfortunately,  this  differential  equation  does  not  have 
a  general  analytic  solution  even  if  11  is  expressed  as  a  simple  dipole  field.  It  is, 
therefore,  common  practice  nowadays  to  study  the  propagation  of  energetic  charged 
particles  in  the  earth  s  magnetosphere  bv  numerical  calculations  using  a  mathe¬ 
matical  model  of  the  earth's  magnetic  field. 

The  numerical  methods  for  the  integration  of  Kq.  (1)  are  well  established  and 

corresponding  computer  programs  have  been  published  by  several  authors 
0  !  0 

(McCracken  ct  al;'  Shea  el  ai  ).  As  is  done  usually  in  cosmic  ray  physics, 
the-  particles  are  characterized  in  these  programs  by  their  rigidity 

it.  Rossi,  R.  ,  and  (Albert,  S.  (1P70)  in  Introduction  to  the  Physics  of  Space, 

MeGraw -Hill,  New  York. 

o.  McCracken,  K.G.,  Rao,  I.R.  ,  and  Shea,  Al.  A.  (IDG  2)  The  Trajectories  of 
t  osmic  Ra vs  in  a  High  Degree  Simulation  of  the  Geomagnetic  Field, 

M.I.T.  lech.  Rpt.  No.  77,  NN  0-2G70. 

ID.  mi,-  ,,  M.  \.  ,  Smart,  D.  1  .  ,  and  Carmichael,  II.  (1D7G)  Summary  of  Cutoff 
Rigidities  Calculated  With  the  International  Geomagnetic  Reference  l,'iel’3 
T 1 1 ■  \  a  rious  Dpi >c bs,  LJKI’  No.  a G  1 ,  A 1-  G  I-~  1  R  -  i  G - 0  I  15,  A 1  )A  028117  8. 
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(2) 
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wnere  c  is  tht'  velocity  of  light.  The  procedures  then  utilize  the  fact  that  Eq.  (1) 
remains  unchanged  if  the  signs  of  q  and  ds  are  reversed.  The  trajectory  of  a 
particle  with  rigidity  1‘  arriving  at  a  specific  location  from  a  specific  direction  is, 
therefore,  determined  backwards  by  tracing  the  path  of  an  identical  particle,  but 
with  opposite  charge,  leaving  that  particular  location  in  the  specified  direction. 

The  computer  simulation  of  particle  trajectories  in  space,  that  is,  the 
t  im  lector-.  -‘  racing  technique,  has  been  used  as  a  basis  for  the  procedure  followed 
in  the  analysis  of  the  effects  of  asymmetric  magnetospheric  current  systems  on 
cosmic  radiation.  According  to  this  procedure,  which  is  illustrated  schematically 
in  figure  1,  a  model  of  the  quiescent  geomagnetic  field  had  to  be  specified  first. 

1  ne  mam  field  near  the  earth  is  normally  represented  by  a  sum  of  spherical  har¬ 
monies  with  coefficients  determined  so  as  to  produce  the  best  fit  to  experimental 
data.  In  reci'ii*  years  great  efforts  have  been  made  to  develop  quantitative  models 
able  to  describe  the  magnetic  field  within  the  entire  magnetosphere.  Tor  the  quiet 
time  magnetic  field  configuration,  several  models  now  exist  which  include  the  con¬ 
tributions  l  rom  magnetopause,  ring  and  tail  currents  in  addition  to  the  main  field 
(Walker1  In  the  present  study,  allowance  had  to  be  made  for  the  fact  that  in  the 
process  of  calculating  a  cosmic  ray  tra.jcctorv  the  magnetic  field  lias  to  be  evaluated 
d  mam  positions  md  that  the  respective  amount  of  computer  time  should  stay 
w  it  nm  reasonable  'units.  1  herefore,  we  hav  e  not  utilized  the  most  .  ophisticated 
models  m  this  studv.  1  or  a  number  of  initial  calculations,  which  were  essentially 
intended  to  give  a  more  qualitative  than  quantitative  overview  >1  the  possible  effects 
ol  local  magnetic  perturbations  on  cutoff  rigidities,  a  simpl"  dipole  representation 
ii  m  been  used  to  speed  up  *  he  computations.  In  ■  billowing.  more  quantitative  analy¬ 
sis,  ’ne  quiescent  m  igncMe  field  was  described  In  the  International  (ieomagnetic 
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effects  of  the  magnetic  perturbations  due  to  the  model  current  systems  have  then 
been  evaluated  by  comparing  the  results  of  trajectory  calculations  using  both 
quiescent  and  perturbed  magnetic  field  models. 


Figure  1.  Schematic  Representation 
of  the  Procedure  Followed  in  the 
Study  of  the  Fffects  of  Asymmetric 
IVIagnetospherie  Current  Systems  on 
Cosmic  Radiation 


Based  on  the  requirements  of  man',  cosmic  ray  studies,  two  quantities  were 
of  special  interest  with  respect  to  the  effects  of  disturbances  in  the  geomagnetic 
field:  the  cutoff  rigidities  of  cosmic  ray  stations  and  the  asymptotic  directions  of 
cosmic  ray  particles. 

The  cutoff  rigidit.  of  ,  specific  location  on  the  earth  and  of  a  specific  direction 
of  incidence  is  generally  defined  as  the  rigidity  below  which  cosmic  ray  particles 
are  inaccessible  ti.  that  lor.ifjon  from  the  specified  direction.  Following  the  pro¬ 
cedure  described  by  Shea  e«  ,1  *  *  cosmic  ray  cutoff  rigidities  have  been  obtained 
b.  determining  for  t  he  cn*  i  re  rigidit-.  spectrum  whether  an  individual  rigidity  had 
■  i  If.  lector,  ii  .  ,  ..  -it)!-  Iron,  infinit-.  diet  is,  had  an  allowed  trajectory  I  or  not .  Fo 
i  given  Im  ‘ion  ■.  Pli  gi  om  igm  t  i(  coordinates  4>  and  lor  a  given  direction, 
ch  i  racten/i-d  t>  tin  /<  nitfi  ingh  U  and  the  azimuthal  angle  o,  (see  Figure  2  for  the 
corresponding  • ;  am  -  of  ivleiviii  ei  calculations  were  initiated  at  a  rigidity  well 
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PC(A,  4>,  a.  O)  the  Lowest  rigidity  for  which  the  trajectory 

calculations  yielded  an  allowed  orbit.  For  the 
past  15  years  this  trajectory -derived  value  has 
been  referred  to  as  the  Stormer  cutoff,  although 
its  calculated  value  is  in  most  cases  well  above 

the  actual  cutoff  rigidity  corresponding  to  the 

1 6 

Stormer  cone  as  defined  by  Stormer; 

<£,  0,  o)  the  effective  cutoff,  given  bv 


P„(A.  $.  0,  o)  =  PK](f  ,  <J>.  0.  o) 


P  (',<!>.  ‘K  o) 
M 


X 


C(P; 


A 


<1>.  0.  0)  dp 


P  (A,  <J>,  0,  ••’> 


whe  re 


C(P;  *.,  0>.  0,  <?) 


1  if  the  orbit  corresponding  to  P  is  allowed 
0  if  the  orbit  corresponding  to  P  is  forbidden. 


As  illustrated  in  Figure  4,  the  asvmptotic  direction  of  a  cosmic  ray  particle 
having  a  specifie  rigiditv  and  arriving  at  a  specific  location  on  the  earth  from  a 
specific  direction  of  incidence  is  the  direction  of  motion  which  this  particle  had  in 
interplanetary  space  prior  to  its  interaction  with  the  earth’s  magnetic  field.  It  is 
obvious  that  asvmptotic  directions  are  a  primary  means  of  relating  cosmic  rav  in¬ 
tensity  variations  at  the  earth  with  the  direction  of  the  particles  in  interplanetary 

p  ]7 

space  (McCracken  et  al  ’  ).  The  t rajeetorv -t  racing  method  enables  an  accurate 

determination  of  asvmptotic  directions  for  a  specific  model  of  the  earth's  magnetic 
1 8 

field  (Shea  and  Smart  ).  In  the  relatively  simple  field  models  used  in  this  analysis 
the  asymptotic  direction  of  a  particle  has  been  determine.!  bv  its  direction  of  motion, 
expressed  in  terms  of  the  geocentric  coordinate  system,  at  a  earth  radii,  where 
the  effects  of  the  geomagnetic  field  on  the  orbit  become  insignificant. 

16.  Stormer,  C.  (1''55)  The  Polar  Aurora.  Oxford  1  imersitv  Press, 

Fair  I.awn,  New  .lersev. 

17.  McCracken.  K.(i.,  Kao,  1  .  K.  .  fouler.  It.  C.,  Shea,  M.  \.  ,  and  Smart.  I  >.!■'. 
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In  order  to  understand  spei  ifii  aspects  of  the  effects  of  asymmetric  mag- 
nctosphcric  current  systems  on  cutoff  rigidities  and  asymptotic  directions  it  was 
also  necessary  to  study  the  detailed  elements  of  particular,  representative  particle 
'  >' a  ici  dories.  l  or  this  purpose  several  computer  programs  were  developed  to  give 
n  a  pin  i  d  renresentat  ions  of  the  partiele  orbits  as  well  as  of  a  number  of  eharac- 
tertstie  parameters  evaluated  along  the  trajectories. 


ASYMPTOTIC 

VELOCITY 


STATION  LONGITUDE  4> 


I  igtire  I.  Iliust  rat  Ion  of  the  I  ad'tttit  ion  of  the  Asvmntotie 
I  i  tv , •:  i '  it:  o i  \  oi i I’na  ch  t  no a>r din j>  '  o  Shea  and  Ste a rt  * '*  I 
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TR.IPI.TH 


An  example  of  a  plot  c  reated  bv  program  TRJPLT2  is  given 
in  Figure  6. 

Program  to  plot  the  trajectory  described  by  a  particle  in 
the  particular  magnetic  meridian  plane  which,  at  any 
'me,  contains  the  particle  and  rotates  with  the  particle 
around  the  geomagnetic  axis.  A  sample  plot  is  given  in 
Figure  7.  This  representation  is  useful  to  study  the 
complexity  of  the  trajectories  and  to  classify  the  orbits 
according  to  their  different  type. 


1'  iglire  a.  Side  and  Top  \  iew,  in 
Geomagnetic  Coordinates,  of  the 
I  ra jeetorv  W  it  hi  n  Geocentric 
Distances  r  '  Sr,,  Described  hv  a 
<  Osm  ie  Hav  Partic  le  \V  ith 
Kigiditv  P  -  1.42  G\  ,  Arriving 
I  rom  the  Vertical  Direction  at 
4  ting  I  ra  n joeh  (geomagnet  i  c 
latitude  -  ij  -  47. 7H5,  geomagnetic 
longitude  <tu  -  ''0.  07 ° ,  as  Calculated 
I  tilr/ing  the  International  Geomag¬ 
netic  Reterenee  Field  Vpproprjate 
'or  Fpoeh  li'tia.  (1 
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SIMPLE  STUDY  ON  THE  EFFECTS  OF  LOCAL  DISTURBANCES  OF 
THE  GEOMAGNETIC  FIELD  ON  COSMIC  RAY  CUTOFF  RIGIDITIES 


In  nnler  to  got  an  insight  into  the  general  aspects  of  the  effects  of  asymmetric 

magnetospheric  current  systems  on  cosmic  radiation,  an  initial  study  was  con- 

luctod  correlating  the  variations  of  the  cutoff  rigidities  of  two  European  neutron 

monitor  stations,  Kiel  and  Jungfraujoch  (see  Table  1  for  corresponding  station  data) 

!h  specific  local  perturbations  of  the  geomagnetic  field.  In  this  analysis  the 

quiescent  geomagnetic  field  was  represented  by  the  field  of  a  simple  dipole  with 

25  2 

magnetic  moment  M  8.06  •  10  ti -cm  ,  having  its  origin  at  the  center  of  the 
earth  and  the  dipole  aKts  intersecting  the  surface  of  the  earth  at  78.  57  ’  N  and 
(io.Tti’W.  i  he  components  of  B  at  geocentric  distance  r,  geomagnetic  latitude  A 
and  geomagnetic  longitude  <|>  were,  therefore,  given  by 


B  t  r.  ,  4>)  -  (2IV  /  r  )  sin  A 

r 


B.tr,  A,  <t) 


(M/r'l  cos  A 


B.  (  r.  A,  <l>)  0 
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In  the  first  part  of  the  studv  the  perturbed  geomagnetic  field  was  modeled  bv 

—* 

superposing  the  disturbance  field.  B1, 


IT l  r. 


4» 


r 


r 


B 


(-4) 


on  the  dipole  field.  Bv  varying  the  adjustable  boundary,  r^,  successively  from 
r  1  eartli  radius  (r  t  to  larger  values  it  was  possible  to  analyze  the  dependence 
of  the  cutoff  changes  on  the  radial  structure  of  geomagnetic  perturbations.  In  order 
to  facilitate  comparisons  of  the  results  with  the  theoretical  work  bv  Treiman’*  on 
the  effect  of  the  equatorial  ring  current  on  cosmic  rav  intensity,  B"  was  chosen  so 
as  to  describe  the  magnetic  field  of  a  current  sheet  flowing  on  a  geocentric  sphere 
with  radius  a  and  having  a  current  densitv  proportional  to  the  cosine  of  the  geomag- 
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The  results  of  these  calculations  are  illustrated  in  Figures  8  and  9.  In 
Figure  8  the  relative  changes  in  the  vertical  main  cutoff  and  in  the  effective  vertical 
cutoff  rigidity  at  Jungfraujorh  are  plotted  as  a  function  of  the  radial  extent,  r^,  of 
a  magnetic  perturbation  according  to  Fq.  (4)  with  B"z(r  £  a)  =  -  100  nT  and  with 
the  values  a  .'i,  4  and  5r(>  for  the  radius  of  the  corresponding  ring  current.  It 
can  be  seen  that  for  r^  5  a  (if  a  =  3r  ),  and  r^  S  4^e  (if  a  -  4re>,  the  dependence  of 
the  relative  cutoff  variation  AP/P  on  r^  can  well  be  approximated  by  a  linear 
function.  Since  AP/P  also  depends  linearly  on  B "  ,  we  find  in  this  range  of  r^ 


BOUNDARY  RADIUS  OF  MAGNETIC  PERTURBATION,  r,  [  re  ] 


Figure  8.  Relative  Changes  in  the  Vertical 
Main  Cutoff  (O),  and  in  the  Kffective 
Vertical  Cutoff  Rigidity  (•)  at  .lunglraujoch 
as  a  Function  of  the  Boundary  Radius,  t'jj, 
of  a  Magnetic  Disturbance  According  to 
Fqs.  (4)  and  (5)  With  B"z(r  '  a)  -100  nT 
and  a  3,  4  atid  5r(,  Superposed  Upon  the 
Main  Dipole  Field  (see  text  for  details).  The 
corresponding  undisturbed  cutoff  rigidities  are 
lJA1  3.  43  CIV,  F  3.  22  CIV 


I  nil  he  rmore ,  the  relative  cutoff  changes  are  almost  independent  of  fj  for 

r  ■  a  lit  a  3r  )  and  for  r  ■  4r  (if  a  '  4r  )  w;th  approximate!  v  equal  amplitudes 

or  ill  values  a  '  Ir  .  This  last  result  is  in  good  qualitative  agreement  with  the 

on.  lusions  obtained  hv  I  reiman  for  the  svmmetric  ring  current.  In  Figure  a 

"iTes|inii ding  results  for  Kiel  obtained  with  the  same  disturbance  Held 

R  It  f  at  lilt)  n  i  and  with  radii  of  the  ririL'  current  equal  to  a  4,  a,  li  and  7  r 
/  e 

ire  shown.  All  hough  these  results  are  in  general  similar  to  those  found  for 
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lungfrau  ioeh,  signi fi i  ant  iiil'lcrcnics  can  he  observed.  First  of  all,  the  dependence 
■)!'  Mu-  relative  cut  off  variation,  _il’/lJ,  .>n  the  boundary  value  r  can  no  longer  be 
approximated  iiv  a  single  linear  function  in  the  ranges  r^  <  a  (if  r^  -  4,  h,  and  6rp) 
and  r.,  ■  fir  iif  a  7  r  >.  On  the  contrary,  we  find 

i>  i  e 


(8) 


<!>) 


for  2.  5r  <  r,,  <  a  (if  a  -  4,  or  fir  )  and  for  2.  T>r  <  fir  (if  a  =  7 r  ).  1  his 

e  15  e  v  e  e 

indicates  that  the  cutoff  rigidity  at  Kiel  is  approximately  five  times  more  sensitive 

p>  perturbations  of  the  /  -component  of  the  geomagnetic  field  at  geocentric  distances 

2.  5r  <  r<fir  than  to  equal  magnetic  disturbances  witnin  1  r  <  r  <  2.  5r  .  It  can  then 
re  e  e 

also  be  seen  from  the  curves  corresponding  'o  a  -  6r(  anti  a  =  7  r  that  the  cutoff 
rmiditv  at  Kiel  is  influenced  bv  changes  In  the  geomagnetic  field  within  geocentric 
li. stances  of  up  to  r  s  8r  in  contrast  to  the  corresponding  value  of  r  *  4r 
found  for  .lungfrau.  |och. 

In  the  second  part  of  the  analysis  the  disturbance  lield  superposed  upon  the 
L'comagncl  ir  dipole  field  was  chosen  to  he 

|l*5"(r,  ' ,  0)  for  0^  *  0  '  0j. 

(10) 

0  elsewhere 

w  hi’  i  i  It  was  again  d'dined  h>  lap  ( a )  a  nd  w  lie  re  4^  and  0^ .  desc  not'  the  geomagnetic 

lon.otii  |c  of  the  western  ni  l  eastern  houndarv,  respectively,  of  the  longitude  sector 

out  tming  the  magnetic  pert nrhation.  Hv  varying  the  parameters  0^.  and  0^,  with 

•  >•  - 1 >i • ,  t  i  s pc .  i  lie  location  on  the  earth  it  was  possible  to  study  the  dependence  of 

die  i  iiio'i  .  hanjes  on  the  longitudinal  structure  of  geomagnetic  perturbations.  Typical 

•'  mils  of  those  calculations  arc  presented  in  figures  10  through  13.  Figure  10 

-how-.  'Ik  relative  .  Iiangi-s  in  the  vertical  main  cutoff  and  in  the  effective  vertical 

utorf  rigiiitv  ii  .!  ung  frarnoi  h  as  a  function  of  the  longitudinal  position  of  the  center 

>•  eridl.iii,  I  <0  ■  0  1  2.  •!  a  magnetic  perturbation  with  constant  width 

0, .  0,.  Ml’,  !!u\  density  |5  ( r  •  ,ti  -  -100  n  I  and  a  -  Ir  .  It  is  evident  that  both 

F  VV  z  e 
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Figure  10.  Relative  Changes  in  the  Vertical  Main 
Cutoff  t<C>>.  and  in  the  Effective  Vertical  Cutoff 
Rigidity  (#)  at  .Jungfraujoeh  as  a  Function  of  the 
Longitudinal  Position,  T  =  +  4^)/2  of  a  Mag¬ 

netic  Disturbance  According  to  Eqs.  (10)  and  (5) 

With  -  90°,  B"z(r  -  a)  =  -100  nT  and 

a  -  4rc,  Superposed  Upon  the  Main  Dipole  Field 
(see  text  for  details).  The  corresponding  undisturbed 
cutoff  rigidities  are;  P^j  =  3.  43  GY,  P  .  -  3.  22  GV 
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Figure  11.  Relative  Changes  in  the  Vertical  Main 
Cutoff  (<}!,  and  in  the  Effective  Vertical  Cutoff 
Rigidity  (•>  at  Kiel  as  a  Function  of  the  Longitudinal 
Position,  T  "  ■  4>|..)  '2  of  a  Magnetic  I  listurhanc  ■ 

According  to  F.qs.  (10T  and  (a)  With  -  lW  90  , 
H"./(c'  a)  -100  nT  and  a  =  4c,,,  Superposed  Upon  the 
Main  Dipole  Field  (see  text  for  details).  The  corres¬ 
ponding  undisturbed  cutoff  rigidities  are;  P  =  1.  fid  G\  , 
P  =  1.T0  G\  ^ 


L pi  to,-  .  If,  ,  tne  '.erti«  1  cutoff  rigiditv  it  .lungfr ..u.ioclv  its-  shov.n  as  a  function  ot 
' ;  a  ■  gnu  t  -  ,gm  t  i  c  long  it  tide  of  1 1 . stern  boundary,  <1>|  ,  ol  the  longitude  sei  tor 
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containing  the  disturbance  field,  with  the  station  being  located  at  the  western 
boundary,  that  is,  with  The  magnetic  perturbation  is  again  characterized 

by  B"z(r  <  a)  =  -100  nT  and  a  =  4rg.  The  data  show  an  extremely  strong  dependence 
of  AP/P  on  the  longitudinal  width  of  the  perturbed  sector,  W  =  -  4>jj,  in  the  range 

0°  <  W  <  30°.  For  W  ~  60°  the  resulting  changes  in  the  cutoff  rigidities  already 
amount  to  ~80%  of  the  final  value  obtained  for  the  W  >  90°.  It  can  be  concluded 
therefore,  that  cutoff  changes  at  midlatitude  stations  are  mainly  determined  by  the 
magnetic  configuration  within  a  longitude  sector  extending  from  the  longitude  of  the 
station  to  approximately  60°  east  of  it.  Figure  13,  finally,  shows  the  corresponding 
results  obtained  for  Kiel.  It  is  again  found  that  both  the  vertical  main  cutoff  as  well 
as  the  effective  vertical  cutoff  rigidity  are  extremely  sensitive  to  magnetic  perturba¬ 
tions  above  and  up  to  -30°  to  the  east  of  the  station.  For  magnetic  disturbances 
with  longitudinal  width  30°  <  VV  <  60°  the  maximum  changes  in  the  cutoff  rigidities 
are  of  the  order  of  -20%.  It  is  interesting  to  note  that  these  changes  are  by  a  factor 
of  ~1.  25  larger  than  the  depressions  resulting  for  \Y  >  90°,  that  is,  for  magnetic 
perturbations  with  a  large  longitudinal  extent. 
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I  igu  re  13.  I  {dative  Changes  in  tin-  \  -  C  ie  il 
Main  Cutoff  (0  ),  and  in  the  Kffectiw  Vertical 
Cutoff  KiyidiK  (•>  at  .1  ungf ra u jocli  .  s  i 
l  'in.  turn  of  the  Ccomagnetir  l  ongitude  of  the 
Paste  rn  Hound  ary,  <I'|  •  of  1  -Magnetic  I'i.-liirh- 
itice  \e.  online  to  laps.  (  10'  end  (a)  With 
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Figure  13.  Relative  Changes  in  the  Vertical  Alain 
Cutoff  (O).  and  in  the  Effective  Vertical  Cutoff 
Rigidity  (•)  ut  Kiel  as  a  Function  of  the  Geomagnetic 
I.ongitude  of  the  Eastern  Boundary,  4>jr,  of  a 
Magnetic  Disturbance  According  to  Eqs.  (10)  and 
(a)  With  <i>\Y  =  B"z(r  '  a)  =  -  100  nT  and  a  =  4re, 
Superposed  Upon  the  Alain  Dipole  Field  (see  text 
for  details).  The  cor- responding  undisturbed  cutoff 
rigidities  are:  P^(  =  1.113  GV,  P  .  =  1.70  CiV 


Based  on  the  results  of  the  simple  study  presented  in  this  section,  the  follow¬ 
ing  main  conclusions  could  be  drawn: 

-  at  low  and  inidlatitude  stations  ( |  A |  <  50°)  changes  in  the  cutoff  rigidities 

arc  mainly  determined  by  variations  of  the  /-component  of  the  magnetic 

field  within  a  longitude  sector  extending  from  the  longitude  of  the  station 

to  approximately  (iO  1  east  of  it  and  within  geocentric  distances 
2  ■ 

r  •  1  r  •  cos  and 
c 

-  .t  locations  with  geomagnetic  latitudes  50°  <  |  A|  <_65°  changes  in 

cutoff  rigidities  are  mostly  due  to  variations  of  the  /-component  of 

the  magnetic  field  within  a  longitude  sector  extending  from  the 

longitude  of  that  particular  location  to  approximately  00°  east  of  it 

■> 

arid  within  geocentric  distances  r  <  10c  •  cos  “a. 


TIIE  EFFECT  OK  A  PARTIAL  KINO  CURRENT  AND  OF  FIELD-ALIGNED  CURRENTS 
ON  COSMIC  RAY  CUTOFF  RIGIDITIES  AND  ON  ASYMPTOTIC  DIRECTIONS 


III  ->  scamd  study,  the  effect  of  a  partial  ring  current  and  tne  effect  of  field - 

aligned  currents  on  the  cutoff  rigidities  at  Kiel,  Jungfraujocn  and  Home  (see  Table  1 

for  corresponding  station  data),  and  on  the  asymptotic  directions  of  cosmic  ray- 

particles  arriving  at  t  icse  locations  w,.s  investigated.  This  analysis  was  based  on 

t  ie  model  current  svstems  suown  in  f  igures  14  and  15,  corresponding  to  the 

dO 

Hostrom  .No.  1  and  2  Models  (Bostrom-  ).  The'  magnetic  fields  of  the  current 

systems  have  been  calculated  according  to  tne  method  proposed  bv  Kisabeth  and 
21 

Kostoker.  Contributions  b\  currents  induced  in  the  ground  were  taken  into  account 
by  approximating  the  earth  as  an  infinitelv  conducting  sphere  covered  by  a  perfect 
insulator  of  depth  D  (D  =  250  km).  The  International  Geomagnetic  Deference  Field 
appropriate  lor  Epoch  1065.0,  in  geomagnetic  coordinates  (Mead^),  was  used  to 
describe  the  quiescent  Ueid;  the  perturbed  model  was  represented  by  the  super¬ 
position  of  the  magnetic  fields  generated  by  the  respective  model  current  systems 
and  the  quiescent  lield. 


I  igure  14.  Side  and  Top  View 
of  Model  Current  System  No.  1 
(partial  ring  current)  and 
Definition  of  its  Descriptive 

I ’a  r  .motors 
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Figure  15.  Side  and  Top  View  of 
Model  Current  System  No.  2  and 
Definition  of  its  Descriptive 
Parameters 


We  first  discuss  the  cutoff  variations  caused  by  a  Model  No.  1  current  system 

g 

with  U  =  90 °,  1.  =  4  r  ,  I  =  2-10  A,  and  an  altitude  of  the  ionospheric  currents 
l(  -  100  km.  According  to  Fukushima  and  Kamide^  such  a  current  system  is  repre¬ 
sentative  of  the  magnetic  disturbances  observed  during  magnetic  storms  in  the 
evening  sector.  It  produces  a  maximum  decrease  of  the  equatorial  magnetic  field 
at  r-  =■  Ir  of  1G0  tiT,  with  the  magnetic  flux  density  B^(r,  A  =  0°,  <I>)  in  the  equatorial 
plane  depending  on  geomagnetic  longitude  and  radial  distance  as  shown  in  Figure  16. 
In  Figure  17a  the  relative  variation  of  the  vertical  main  cutoff,  P  ,  and  of  the 
effective  cutoff,  p  ,  as  a  func  tion  of  the  longitude  of  the  center  meridian  of  the 
current  system.  T,  is  shown  for  Kiel.  The  maximum  effect  with  the  main  vertical 
cutoff  lowered  h\  24.  5'A  and  the  effective  vertical  cutoff  decreased  by  26.  5lr  is 
found  for  T  120°,  that  is,  when  the  station  is  located  near  the  western  boundary 
of  the  longitude  region  with  maximum  perturbation  of  the  equatorial  magnetic  field. 

In  Figure  17b  corresponding  results  ere  shown  for  the  mid-latitude  station 
Jungfrau  joch.  It  is  obvious  that  especially  the  variation  of  the  effective  cutoff  shows 
a  much  more  pronounced  dependence  on  T  at  Kiel  in  comparison  with  .1  ungf raujoch. 
This  result  is  in  agreement  with  the  fact  pointed  out  by  Dorman'5  that  the  typical 
cosmic  ray  increases  observed  during  geomagnetic  ally  active  periods  depend  more 
clearly  on  longitude  for  stations  with  small  geomagnetic  threshold  than  for  low - 
I  ctitude  stations.  The  latitudinal  struc  ture-  of  the  cutoff  changes  caused  by  a 
Model  \o.  1  current  end  its  dependence  on  the  radius  of  the  par4i al  ring  current 

are  lllust  r  .ted  in  gable  2.  Ilc-rc-  the-  relative  variation  of  the  effective-  vertical 
cutoff  rigidit  .  normalized  to  tin  n,  iximum  change  in  tin-  cquatorj  il  surf  ice  mag¬ 
netic  field  mil  iMTip-cl  "Ut  the  iii'i  r-.  il  JO  '  T  ’  120  ,  is  listed,  lor  'In-  stations 

Kiel,  J  ungf  r  -  a  |oc  li  ,  rid  Koine.  |t  is  e\  ulc  i  it  that  t  Ik-  cutoff  v  in  ih"iis  at  Kiel  depend 
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significantly  on  the  radius  L  whereas  at  Jungfrnujoch  and  Home  they  do  not.  The 
resulting  latitude  effect  is  a  characteristic  property  of  the  particular  geometrical 
structure  of  the  current  system. 
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Figure  16.  Magnetic  Flux  Density  IMr,  -  0l>,  4>' 
in  the  Fquatoeial  Plane  of  a  Current  S'vstem  Accord¬ 
ing  to  Model  No.  1  (partial  ring  current)  with  T  -  0  ", 
\V  -  DO'’,  !.  4r(,  and  i  -  2-  lO^A,  as  a  Function  of 

Geomagnetic  Fongitude  and  Kadial  Distance 


Figure  17.  Longitudinal  LITect  of 

Mod(d  Current  .svsfein  No.  1 

\\  it h  W  DO",  1.  Ir,.  and  I  2-  10  \. 

on  the  \ertieal  Main  Cutoff  iO>  mil 

on  tin  | ; fleet  ive  \  ■  Died  (  ut  oi  l 

I  f  i  £_■  i  d  i  t  \  (  •>  at  Kiel  aid  .1  ungfr  ■  u  joe ! . 
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Table  2.  Effect  of  Model  Current  System  No.  1  (partial  ring  current)  With 
W  =  90°  on  the  Effective  Vertical  Cutoff  Rigidity  at  Kiel,  Jungfraujoch  and 
Rome,  Averaged  for  90°  s  T  £  120°  and  Normalized  to  the  Maximum  Change 
in  the  Equatorial  Surface  Magnetic  Field 


Station 

Effective  Vertical 
Cutoff  Rigidity 

UGRF,  Epoch  1965.  0) 

Normalized  Relative  Change  in 
Effective  Vertical  Cutoff  Rigidity 
[To  /  100  nT] 

[GV] 

L  =  3r 

e 

H 

0) 

L  =  5r 

e 

Kiel 

2.  26 

8.  5  ±  0. 5 

15.  8  ±  0.  5 

19.  8  ±  0.  6 

Jungfraujoch 

4.  47 

4.  6  ±  0.  5 

5.  3  ±  0.  5 

4.  3  ±  0.  6 

Home 

6.  16 

3.5  ±  0.  2 

3.  5  ±  0.  2 

3.4  i  0.  2 

The  effect  of  a  partial  ring  current  on  asymptotic  directions  is  illustrated  in 
Figures  It)  and  19.  Both  figures  refer  to  cosmic  ray  particles  arriving  from  the 
vertical  direction  at  Jungfraujoch,  and  to  a  current  system  with  the  parameters 
f  =  90°,  W  =  90°,  F  =  4t'e  and  1  =  2-  10^A.  In  Figure  1!),  the  asymptotic  directions 
obtained  for  specific  rigidities  in  the  perturbed  geomagnetic  field  model  are  com¬ 
pared  to  corresponding  results  obtained  by  using  the  KIRI-',  Epoch  1965.0.  In 
Figure  19  the  absolute  changes  of  the  asymptotic  directions  in  longitude  and  latitude 
are  plotted  as  a  function  of  the  rigidity  of  the  particles.  It  can  be  seen  that  for 
rigidities  F  >  5.5  GV  the  effect  is  mainly  restricted  to  a  westward  shift  .n  longitude 
with  amplitudes  |Ai//|  <  25°.  As  the  rigidity  decreases  below  F  *  5.5  GV,  the 
magnitude  of  the  effect  increases  strongly,  reaching  maximum  changes  in  both 
asymptotic  longitude  and  latitude  for  rigidities  just  above  the  value  of  the  quiescent 
main  cutoff.  It  is  interesting  to  note,  however,  that  the  asymptotic  direction 
corresponding  to  the  respective  main  cutoff  rigidity  remains  practically  unchanged, 
regardless  of  the  actual  value  of  the  main  cutoff  rigidity. 

Due  to  five  fact  that  tin-  effects  of  magnetospheric  current  systems  on  cutoff 
rigidities  and  on  asymptotic  directions  are  closely  correlated,  all  the  results 
desci'ibing  the  dependence  of  the  cutoff  changes  on  the  parameters  of  the  current 
system  apply  analogously  to  \  a  ri  hions  of  tin*  asymptotic  directions. 
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Figure  18.  Effect  of  a  Current  System  According  to 
Model  No.  1  (partial  ring  current.  PRC)  With  T  =  90° 
W  =  90°,  1.  =  4r(,  and  I  =  2*  lO^A,  on  the  Asymptotic 
Directions  for  Cosmic  Ray  Particles  of  Vertical 
Incidence  at  Jungfraujoch  (+) 


>-  60” 


Figure  19.  Changes  in  Asymptotic 
1  attitude,  A,  and  Longitude,  i//,  of 
Cosmic  Ray  Particles,  Incident  l-  rom 
the  Vertical  Direction  at  Jungfraujoch, 
Due  to  a  Current  System  According  to 
Model  No.  1  (partial  ring  current,  PRC) 
With  T  =  90° ,  W  =  90°,  I.  =  4  r  ,  and 
I  =  2-  106A  l' 


i  lii'  effects  of  field-aligned  currents  on  cosmic  radiation  are  in  general  less 
pronounced  compared  to  the  effects  caused  by  a  partial  ring  current.  The  longi¬ 
tudinal  effect  of  a  model  current  system  according  to  Figure  la  is  illustrated  in 
Figure  20,  fur  Hie  parameters  1.1  -  4r  ,  1.2-  5.  5r  ,  1-  2-10^A  and  II-  100  km. 

e  e 

In  this  plot  the  relative  variation  of  the  main  cutoff  and  of  the  effective  vertical 
cutoff  at  Kiel  are  shown  as  a  function  of  the  longitudinal  position  of  the  current 
system.  I  or  tile  specified  directions  of  the  curren*  vectors,  the  cutoff  rigidities 
are  m,  re  sell  v.  hen  Hie  station  is  located  more  than  -  5  ’  west  of  the  current  system. 
I  or  Hie  m  on  cutoff  i  maximum  increase  of  :i.  occurs  for  T  117°  whereas  the 
effect  r.r  i  utofi  rigidity  shows  a  maximum  increase  of  4 .  i!  at  T  ■■  10a"’.  The  cutoff 

rigidities  iii-  decreased  when  the  station  is  situated  east  of  the  current  system. 

1  he  m  , Minimi  depression  of  the  m  i in  ciiHiff  is  found  for  T  -  !'0J  with  an  amplitude 


of  8.4%.  For  the  effective  cutoff  a  maximum  decrease  of  5.7%  results  at  T  =  82°. 
A  similar,  but  less  pronounced  longitudinal  structure  is  found  at  midlatitudes  indi¬ 
cating  that  the  amplitude  of  the  effect  depends  strongly  on  latitude.  The  latitude 
dependence  of  the  cutoff  changes  caused  by  a  Model  No.  2  current  system  is  illus¬ 
trated  in  Figure  21.  Here  the  relative  variation  of  the  effective  vertical  cutoff 
rigidity  for  Kiel,  Jungfraujoch  and  Home  averaged  over  the  interval  80°  S  T  i  90° 
and  normalized  to  the  current  intensity  is  plotted  as  a  function  of  the  difference 
A  '-  between  the  geomagnetic  latitude  corresponding  to  the  foot  point  of  the  inflowing 
field  aligned  current  [  *  j  =  are  cos  (r  /LI)®’  j  and  the  geomagnetic  latitude  of  the 
stations.  ’.  The  different  data  points  have  been  obtained  for  LI  =  3.0,  4.0,  5.0  r 
and  1.2  =  4.5,  5.5,  6.5  r  respectively.  As  can  be  seen  the  resulting  correlation 
can  be  approximated  by  an  exponential  of  the  form  fexp(-A2/5.  14®)). 
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Figure  20.  Longitudinal  Effect  of 
Model  Current  System  No.  2  With. 
1.1  =  4re,  1.2  -  5.  5re  and  I  =  2-  106A 
on  the  V  ertical  Main  Cutoff  (  O*  and 
on  the  Effective  Vertical  Cutoff 
Rigidity  (•)  at  Kiel 
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Figure  21.  Average  Effect  of 
Model  Current  System  No.  2  on  the 
Effective  Vertical  Cutoff  Rigidity 
at  (ieomagnetic  Longitude  <t>  05° 

for  80°  •  T  '  90° 
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Due  to  the  fact  that  in  this  analysis  the  111  111  lias  been  used  to  model  the 
quiescent  geomagnetic  field,  the  quantitative  results  obtained  can  actually  be  re¬ 
lated  to  real  experimental  measurements  and  observations.  In  particular,  a 
comparison  of  theoretical  and  experiment  i]  data  will  allow  us  to  verify  and  adjust 
the  geometry  and  mtensit.  of  the  main  current  systems  used  to  model  the  disturbed 
magnetospho  re. 
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in  tlu  geomagnetic  field  and  their  effects  on  cutoff  rigidities  and  asymptotic  direc¬ 
tions.  In  particular,  our  calculations  show  that  the  method  which  was  up  to  now 
generally  used  by  cosmic  ray  physicists  to  calculate  changes  in  cutoff  rigidities 

from  the  average  worldwide  depression  of  the  equatorial  surface  magnetic  field 

24 

during  geomagnetic  storms,  that  is,  from  the  parameter  D  ,  (Obayashi; 

25  ST 

Kliickiger  et  al“‘  ),  can  lead  to  substantial  errors  of  the  order  of  up  to  100%.  The 

proper  analysis  lias  to  take  into  account  the  asymmetric  longitudinal  structure  of 

the  magnetic  perturbation.  A  method  which  makes  allowance  for  this  effect  will  be 

proposed  in  another  report. 

Of  further  significance  is  also  the  fact  that  the  present  study  clearly  indicates 
that  in  order  to  be  able  to  describe  the  propagation  of  cosmic  ray  particles  of 
rigidities  1’  <  2  Cl\  through  the  domain  of  the  earth's  magnetic  field  in  full  con¬ 
sistence  with  the  observations,  even  the  model  of  *he  quiescent  geomagnetic  field 
should  include  the  contributions  due  to  external  sources.  The  technique  presented 
in  this  report  allows  for  this  effect  and  can,  therefore,  be  essentially  useful  in  the 
study  of  the  origin  of  particles  incident  at  tile  earth  at  high  latitudes,  especially  in 
the  vicinity  of  open  field  lines.  This  problem  is  still  not  solved  satisfactorily 

although  it  is  of  considerable  interest  for  nianv  applications  (IJcbrunner  and 
2(i 

Ian  kuoo.l  |_  f  inally,  from  an  inspection  of  all  the  results  obtained  in  this  study 
re  I  ding  the  changes  in  cutoff  rigidities  and  in  asymptotic  directions  with  specific 
m  ignefic  perturbations,  and  from  the  comparison  of  the  magnitude  of  these  dis¬ 
turb.  1 1 it -i  s  i s  a  function  of  geocentric  distance,  geomagnetic  latitude  and  longitude 
1 1 1 1  the  actual  magnitude  of  the  external  contributions  to  the  geomagnetic  field  it 
c  >n  be  concluded  that  from  the  point  of  view  of  cosmic  ray  research  for  the  re  pro - 
■  illation  of  the  main  geomagnetic  field  ‘here  is  probably  no  need  to  include  spherical 
a.  rn  i>  mi.  a-,  with  orders  much  higher  than  nine. 

In  ■  ainm  iri/ing,  the  new  insight  obtained  in  this  study  into  the  effect  of  perturba¬ 
tion  'it  'he  geimi  ignetii  field  mi  the  propagation  of  cosmic  ray  particles  near  the 
>. '  ti'tli  .'.ill  a  rt.ainl  i  liable  m  impro'.  ed  .nilvsis  of  cosmic  ray  data  obtained  during 
magnetic  ill1,  •••five  periods  and  hopefully  lead  to  an  increased  use  of  cosmic  ray 
inoi.-uri  men*.-'  is  m  .dditioiial  1 1 ,< >  1  in  the  study  of  the  disturbed  magnetosphere. 

2  1.  I  >b..y  .sin.  T.  ( 1  **"»;»>  1  hit ry  of  high  energy  particles  into  the  sola  r  ionosphere, 
li'pt.  Ionosphere  and  Space  lies,  .lupati,  l  .'KXo.  2):  201. 

- - - |/VV\ 

2a.  I  luckiger.  11.,  I  >■  hrnmn •  r,  II.,  \  reus,  M.,  md  Hinder,  (>.  (1275)  Cutoff 

rigidity  i  nation.-,  of  l.urope  >n  mid- latitude  stations  during  the  September  1074 

I  orhnsii  dec  re  i  a  ,  i  •  >nfe  reiiee  Papers,  IdthlCRC,  .Munich,  4:1221. 

- : -  vAA^. 

2n.  I  ah  rutuie  r.  II.,  >nd  I  ..  >.  k.M  mil,  .1.  A.  Ilh.'IO)  The  spatial  anisotropy,  rigidity 

spectrum,  ,  nd  p  rop  :g  it  n  >n  cii.  i  r  a  ct  e  rist  ies  of  the  relativistic  sol.  >  r  pa  rt  ides 

during  tile. "vent  m  M.ev  V,  I'oil,  .1.  <  ieoplivs.  lies.  li.u :  (>  ha  2. 
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